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ABSTRACT 

We have constructed the first all-sky CMB temperature and polarization lensed maps 
based on a high-resolution cosmological A^-body simulation, the Millennium Simula- 
^ ' tion (MS). We have exploited the lensing potential map obtained using a map-making 

procedure (Carbone et al. 2008) which integrates along the line-of-sight the MS dark 
matter distribution by stacking and randomizing the simulation boxes up to z = 127, 
' and which semi-analytically supplies the large-scale power in the angular lensing po- 

tential that is not correctly sampled by the iV-body simulation. The lensed sky has 
been obtained by properly modifying the latest version of the LcnsPix code (Lewis 
2005) to account for the MS structures. We have also produced all-sky lensed maps 
of the so-called V'-E and ipB potentials, which are directly related to the electric and 
magnetic types of polarization. The angular power spectra of the simulated lensed tem- 
perature and polarization maps agree well with semi- analytic estimates up to ^ < 2500, 
while on smaller scales we find a slight excess of power which we interpret as being due 
to non-linear clustering in the MS. We also observe how non- linear lensing power in 
the polarised CMB is transferred to large angular scales by suitably misaligned modes 
in the CMB and the lensing potential. This work is relevant in view of the future 
I CMB probes, as a way to analyse the lensed sky and disentangle the contribution 

. from primordial gravitational waves. 
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1 INTRODUCTION 

The Cosmic Microwave Background (CMB) is character- 
ized both by primary anisotropics, imprinted at the last 
scattering surface at redshift z ^ 1100, and by secondary 
anisotropies caused along the way to us by density inhomo- 
geneities and re-scattering off electrons that are freed during 
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the epoch of reionization, and heated to high temperature 
when massive structures virialize. 

On one hand, the primary CMB anisotropies give direct 
insight into the structure of the very early Universe, and are 
one of the principal pillars on which the standard cosmolog- 
ical ACDM model is founded. The temperature anisotropy 
power spectrum has now been measured to very high preci- 
sion (e.g. Komatsu et al. 2008) yielding tight constraints on 
the basic parameters of the cosmological model. 

On the other hand, one of the most important mecha- 
nisms that can generate secondary anisotropies is the weak 
gravitational lensing of the CMB, which arises from the dis- 
tortions induced in the geodesies of CMB photons by gra- 
dients in the gravitational matter potential (Bartelmann & 
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Schneider 2001; Lewis & Challinor 2006). The remapping of 
points produced by lensing induces non-Gaussianities in the 
observed CMB sky, and also changes the power spectra of 
the perturbations. 

The CMB is also expected to be polarized at the ~ 10% 
level principally because of Thomson scattering of photons 

05 free electrons during recombination. Thomson scattering 
generates hncar polarization only, which is usually expressed 
in terms of the Stokes parameters Q and U. They can in 
turn be decomposed into coordinate independent E- and 
S-modes of polarization (e.g. Zaldarriaga & Seljak 1997) 
with opposite parities, the so-called "electric" and "mag- 
netic" typos of polarization. To linear order in perturbation 
theory, primordial scalar (density) perturbations can only 
generate S-polarization, while primordial vector and tensor 
(gravitational waves) perturbations can generate both scalar 
E- and pseudoscalar B-polarization (Seljak & Zaldarriaga 
1997, e.g.). In particular, primary B-modc polarization rep- 
resents the imprints left from the primordial gravitational 
waves (GWs) on the CMB (Hu et al. 1998; Kamionkowski 

6 Kosowsky 1999; Hu & Dodclson 2002): if initial fluctu- 
ations are created very early, e.g. during inflation so that 
the vector growth is damped, primary B-modes are pro- 
duced only by tensor perturbations that, being damped at 
last scattering by the horizon entering, produce the largest 
amount of temperature quadrupolc anisotropy and, conse- 
quently, by Thomson scattering, the largest amount of polar- 
izatiou (Pritcliard & Kamionkowski 2004). The primordial 
gravitational radiation is thought to have been generated 
by quantum fluctuations of the metric tensor during the in- 
flationary era, with a strain amplitude proportional to the 
square of the inflation energy scale. Consequently, the in- 
direct detection of this relic gravitational background via 
the direct observation of the primary B-mode polarization, 
as expected from future dedicated CMB missions by ESA 
and NASA^ , will shed light on the physics of the very early 
Universe and will represent a powerful magnifier on the in- 
flationary era and the very first moments in the existence 
of the Universe (Kamionkowski et al. 1997; Zaldarriaga & 
Seljak 1997; Seljak & Zaldarriaga 1997). 

However, as for the CMB temperature, there are mech- 
anisms also for the polarization that can produce secondary 
B-modes, with the dominant one being again gravitational 
lensing, i.e. cosmic shear (CS), which distorts the pri- 
mary CMB pattern, in particular converting E- into B- 
polarization (Zaldarriaga & Seljak 1998), even in case of ab- 
sence of primary B-modes. Although comparable, B-modes 
from primordial GWs exhibit their peak at multipoles I « 
100, corresponding to the degree scale, while, for lensed B- 
modes, the peak is at / ~ 1000, corresponding to the ar- 
cminute scale. Nonetheless, if the energy scale of infiation 
is 1/^/* < 4 X lO^'^GeV, the CS-induced curl represents a 
foreground for the i « 50 — 100 primordial GW-induced pri- 
mary B-polarization (Cabella & Kamionkowski 2005). This 
could limit the extraction of the gravity wave signal if not 
taken into account correctly (e.g. Seljak & Hirata 2004), 
even though forthcoming CMB probes will have in princi- 
ple the sensitivity and the instrumental performance for the 
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detection of the CMB anisotropies in total intensity and po- 
larization. 

A precise knowledge of the lensing effects would also 
provide new insights and constraints on the expansion his- 
tory of the Universe, on the process of cosmological structure 
formation (Acquaviva & Baccigalupi 2006; Hu et al. 2006) 
and on the cosmological parameter estimation (Smith K. et 
al. 2006; Smith S. et al. 2006). In particular, for a correct 
interpretation of the data from the forthcoming Planck satel- 
lite^, it will be absolutely essential to understand and model 
the CMB lensing, as the satellite has the sensitivity for mea- 
suring the CMB lensing with good accuracy. We note that 
a first detection of CMB lensing in data from the Wilkin- 
son Microwave Anisotropy Probe (WMAP"*) combined with 
complementary data has already been claimed by Smith et 
al. (2007) and Hirata et al. (2008), and evidence for weak 
gravitational lensing of the CMB has been observed at >3- 
sigma significance by Reichardt et al. (2008). 

From the arguments above it follows that the next gen- 
eration CMB experiments will require a detailed lensing re- 
construction and an accurate de-lensing methodology. One 
can try to reconstruct the gravitational lensing effects using 
the so-called quadratic and maximum-likelihood estimators 
(Hu & Okamoto 2002; Hirata & Seljak 2003), which allow 
to reproduce at some level of precision the lensing poten- 
tial from the observed CMB itself, and to invert the pho- 
ton geodesic remapping induced by the cosmic shear. Up to 
now these methodologies have been applied only to limited 
patches of sky (Amblard et al. 2004) and/or under the hy- 
pothesis of a Gaussian distribution of the lensing sources 
(e.g. Seljak & Hirata 2004). However, this is only a first 
order approximation, since the non-linear evolution of the 
cosmic structures induces non-Gaussian features in the lens- 
ing potential, features that in turn have an impact on the 
non-Gaussian statistic of the CMB produced by the point 
remapping caused by the cosmic shear itself on the CMB 
fluctuation pattern. 

It is thus very interesting to test the performance of 
these estimators on full-sky, lensed CMB maps which do 
include the effects of the non-linear structure evolution at 
all the orders. This demands detailed simulated lensed CMB 
maps. 

The increasing availability of high-resolution iV-body 

simulations in large periodic volumes makes it possible to 
directly simulate the CMB distortions caused by weak lens- 
ing using realistic cosmological structure formation calcular 
tions. Our previous work (Carbone et al. 2008) represents 
a first step into this direction. Indeed, existing studies al- 
ready give access to statistical properties of the expected 
all-sky CMB lensing signal (see e.g. Lewis 2005, and refer- 
ences therein) , but these studies are based on 'semi-analytic' 
calculations that use approximate parameterizations of the 
non-linear evolution of the matter power spectrum. On the 
other hand, up to now A'^-body numerical simulations have 
been used to lens the CMB only in limited patches of sky 
(Das & Bode 2008), or to produce full-sky convergence maps 
confined to low redshifts a ~ 1 (Fosalba et al. 2007; Teyssier 
et al. 2008). 
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In Carbone et al. (2008) we have developed and de- 
scribed a procedure which gives access to the full statistics 
of the lensed CMB signal, including non-linear and non- 
Gaussian effects on the full-sky. This should allow improve- 
ments in the methods for separating the different contribu- 
tions to CMB anisotropics in the data, which would help 
substantially to uncover all the cosmological information in 
the forthcoming observations. 

In this paper we apply the methodology developed in 
Carbone et al. (2008) to the construction of all-sky, lensed 
simulated temperature and polarization maps. In Section 2 
we describe the procedure utilized to lens, which is based 
on the simulated lensing potential and deflection angle tem- 
plates obtained in Carbone et al. (2008). In Section 3 we de- 
scribe the resulting simulated lensed maps. In Section 4 and 
Section 5 we apply statistical analyses to the obtained tem- 
perature and polarization lensed maps showing consistency 
and possible differences with respect to semi-analytical ex- 
pectations. Finally in Section 6 we draw the conclusions of 
this work and outline next steps and future applications. 



2 CMB LENSING THROUGH THE 
MILLENNIUM SIMULATION 

Weak lensing of the CMB deflects photons coming from an 
original direction n' on the last scattering surface to a direc- 
tion n on the observed sky, so a lensed CMB field is given 
by ^(n) = ^(n') in terms of the unlensed field X = T,Q,U 
(e.g. Lewis 2005). The displacement of the points is deter- 
mined by the integral of the gravitational potential along 
the line of sight to the last scattering surface, as we review 
below. 

In what follows we will consider only the small-angle 
scattering limit, i.e. the case where the change in the co- 
moving separation of CMB light-rays, owing to the deflec- 
tion caused by gravitational lensing from matter inhomo- 
geneities, is small compared to the comoving separation be- 
tween the undeflected rays. In this case it is sufficient to cal- 
culate all the relevant integrated quantities, i.e. the so-called 
lensing-potential and its angular gradient, the deflection- 
angle, along the undeflected rays. This small-angle scattering 
limit corresponds to the Born approximation. 

Adopting conformal time and comoving coordinates in 
a flat geometry (Ma & Bertschinger 1995), the integral for 
the projected lensing-potential due to scalar perturbations 
with no anisotropic stress reads 



10- 



*(n) 



r <l>(r-n; 770 — r) 



dr , 



(1) 



where r is the comoving distance, r* ~ 10* Mpc is its value 
at the last-scattering surface, 770 is the present conformal 
time, and $ is the physical peculiar gravitational potential 
generated by density perturbations (Hu 2000; Bartelmann 
& Schneider 2001; Refregier 2003; Lewis & Challinor 2006). 

Actually, the lensing potential is formally divergent ow- 
ing to the 1/r term near r = 0; nonetheless, this divergence 
affects the lensing potential monopole only, which can be set 
to zero, since it does not contribute to the deffection-angle. 
In this way the remaining multipoles take a flnite value and 
the lensing potential field is well defined (Lewis & Challi- 
nor 2006). The vector n' is obtained from n by moving its 
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Figure 1. Top panel: The dotted orange line represents the simu- 
lated lensing potential power spectrum obtained via line-of-sight 
integration across the MS dark matter distribution up to 2: = 127. 
The blue solid line is the same as the orange dotted one after rein- 
stating the large scale power with the use of the LS-adding tech- 
nique. The dot-dashed black line represents the power spectrum of 
the lensing potential obtained with the CAMB code stopping the 
line-of sight integration ai z = 127, and including an estimate of 
the non-linear contributions (Smith et al. 2003). The long-dashed 
light green line is the same as the dot-dashed black one with 
line-of-sight integration up to z = 1100. Finally, the dashed vi- 
olet line represents the linear lensing potential power spectrum 
from the CAMB code in the linear approximation and integrat- 
ing up to 2 = 127. Bottom panel: The ratios between the power 
spectra shown in the top panel and the CAMB linear lensing po- 
tential spectrum up to 2 = 127. It is worth to note that there is 
no difference using z = 1100 or z = 127 for the semi-analytical 
expectations. 



end on the surface of a unit sphere by a distance |Vn^(n)| 
along a geodesic in the direction of Vn*I'(n), where [l/r]Vn 
is the two dimensional (2D) transverse derivative with re- 
spect to the line-of-sight pointing in the direction fi = (i9, ip) 
(Hu 2000; Challinor & Chon 2002; Lewis 2005). We assume 
Vn4'(n)| to be constant between n and n', consistent with 
the Born approximation. 

If the gravitational potential $ is Gaussian, so is 
the lensing potential. However, the lensed CMB is non- 
Gaussian, as it is a second order cosmological effect produced 
by cosmological perturbations onto CMB anisotropics, yield- 
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MS T difference map 




Figure 2. Top panel: Difference T map between the lensed and unlensed T fields obtained by supplying the LensPix code with the 
spherical harmonic coefficients extracted from the MS lensing potential map and by implementing in it the LS-adding technique, as 
described in the text. Units in mK. A histogram equalized color mapping has been used to increase the contrast. Middle panel: Modulus 
of the polarization AP = ^ AQ^ _|_ At/^^ where AQ and bJJ are the difference Q and U maps obtained using the same technique as for 
the T difference map of the top panel. Units in mK. Bottom panel: Map of the deflection-angle modulus obtained as angular gradient of 
the lensing-potential whose power spectrum is represented by the solid blue line in Fig. 1 
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ing a finite correlation between different scales and thus non- 
Gaussianity. This is expected to be most important on small 
scales, due to the non-linearity already present in the under- 
lying properties of lenses. 

In order to generate full-sky T, Q, U maps lensed by the 
matter distribution of the Millennium Simulation, we have 
modified the publicly available LensPbc code* (LP), which 
is described in Lewis (2005). In its original version, in fact, 
this code lenses the primary CMB intensity and polariza- 
tion fields via a Gaussian realization, in the spherical har- 
monic domain, of the lensing potential power spectrum as 
extracted from the publicly available Code for Anisotropics 
in the Microwave Background (CAMB'"'). Our modification 
(hereafter referred to as "MS-modified-LP") consists in forc- 
ing LP to deflect the CMB photons using the fully non- 
linear and non-Gaussian lensing potential realization ob- 
tained from the MS using the procedure briefly summarized 
here, which was presented by Carbono et al. (2008); we refer 
the reader to that paper for more detail. 

The MS is a high-resolution iV-body simulation for a 
ACDM cosmology consistent with the WMAP flrst year re- 
sults (Spergel et al. 2003), carried out by the Virgo Con- 
sortium (Springel et al. 2005). It uses about 10 billion col- 
lisionless particles with mass 8.6 x 10*/i~^Mo, in a cubic 
region 500/i~^Mpc on a side which evolves from redshift 
z, = 127 to the present, with periodic boundary conditions. 
Our map-making procedure is based on ray-tracing of the 
CMB photons in the Born approximation through the three- 
dimensional field of the MS peculiar gravitational potential. 
In order to produce mock lensing potential maps that cover 
the past light-cone over the full sky, we stack the peculiar 
gravitational potential grids around the observer (which is 
located at z = 0), exploiting the pre-computed and stored 
snapshots of the simulation. The spacing of the time out- 
puts of the MS simulation is such that it corresponds to 
an average distance of 140/i^^Mpc (comoving) on the past 
light-cone. We fully exploit this time resolution which, at 
high accuracy, allows to avoid the adoption of time inter- 
polation techniques, and use all the 63 outputs of the sim- 
ulation along our integration paths. In practice this means 
that the data corresponding to a particular output time is 
utilized in a spherical shell of average thickness 140 h~^Mpc. 
Moreover, the total volume around the observer up to «* is 
divided into spherical shells, each of thickness 500/i~^Mpc. 
All the MS boxes falling into the same shell are translated 
and rotated with the same random vectors generating a ho- 
mogeneous coordinate transformation throughout the shell, 
while the randomization vectors change from shell to shell. 
The peculiar gravitational potential at each point along a 
ray in direction n is spatially interpolated from the pre- 
computed MS grid which possesses a spatial resolution of 
about 195/i~^kpc. The deflection angle is computed along 
the line of sight as well, by numerically evaluating the grav- 
itational potential gradient and interpolating at each point 
along the line of sight (Carbone et al. 2008). 

Being repeated on scales larger than the box size, the re- 
sulting weak lensing distortion lacks large scale power, which 
manifests itself in the lensing potential power spectrum as 
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an evident loss of large scale power with respect to semi- 
analytic expectations, which is most noticeable for multi- 
poles smaller than I ~ 400. This has been cured by aug- 
menting large scale power (LS-adding) directly in the angu- 
lar domain, a procedure which we exploit here as well. More 
speciflcally, we have implemented the LS-adding technique 
directly into the LensPix code as we now explain. We have 
again split the spherical harmonics domain into two multi- 
pole ranges; < I < 400, where the MS fails in reproducing 
the correct lensing potential power due to the limited box- 
size of the simulation, and / > 400, where instead the power 
spectrum is reproduced correctly by the Millennium Simu- 
lation (sec Fig. 1). On the latter interval of multipolcs, we 
have extracted the corresponding ensemble 'iff^ of lensing- 
potential spherical harmonic coefflcients produced by the 
MS lens distribution. We have modified the LensPix code 
so that it reads and uses these MS harmonic coefficients 
on the corresponding range of multipoles. On the interval 
< Z < 400, instead, we let LensPix generate its own en- 
semble of spherical harmonic coefficients which are a 
realization of a Gaussian random field characterised by the 
CAMB semi-analytic lensing-potential power spectrum (in- 
cluding the estimate of the contribution from non-linearity 
Smith et al. 2003) inserted as input in the parameter file of 
LP. 

Since on low multipolcs the effects of the non- 
Gaussianity from the non-linear scales are negligible and 
the ^^im are independent, every time that we run the MS- 
modified-LP, we generate a joined ensemble of ^(m, where 

= ^fLP for < Z < 400 and ^im = for I > 400. 

This technique achieves two goals: firstly it reproduces 
correctly the non-linear and non-Gaussian effects of the MS 
non-linear dark matter distribution on multipoles I > 400, 
including at the same time the contribution from the large 
scales at Z < 400, where the lensing potential follows mostly 
the linear trend as shown from the light-green dot-dashed 
line in Fig. 1. Secondly, it allows to take correctly into ac- 
count the cross-correlation between the temperature and the 
lensing-potential (e.g. Lewis & Challinor 2006) due to the 
Integrated Sachs- Wolfe (ISW) effect on the low multipoles, 
i.e. on the largo scales (this effect is instead negligible at 
Z > 200 Afshordi 2004). In fact, as it happens for the semi- 
analytic contribution to the lensing-potential, also the pri- 
mary temperature and polarization fields are generated by 
LP as Gaussian realizations, in the harmonic domain, of 
the corresponding power spectra obtained by CAMB, which 
correctly includes the contribution to the temperature from 
the ISW effect. We emphasise that we always run CAMB 
using the same cosmological parameters as the MS speci- 
fied in Section 2, and, for consistency with the MS map- 
making procedure, we have fixed the maximum redshift of 
the line-of-sight integration for the CAMB lensing sources 
at Zmax = 127. Indeed, the lensing power from even higher 
redshifts is negligible for CMB lensing, as we show in Fig. 1 
where the light-green long-dashed line overlaps the black 
dashed-dotted line perfectly. 

On multipoles Z > 400, the cross-correlation between 
the lensing-potential and the temperature is negligible, even 
if there could be some residual contribution coming from 
the non-linear Rees-Sciama effect. In this work we do not 
consider the cross-correlation due to this second order effect. 

To generate the lensed T, Q, U fields from the MS- 
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Figure 3. Difference maps between the ipE lensed and unlensed fields, obtained in the MS-modified-LP case (upper panel) and in the 
unmodified-LP case (lower panel). Units in mK. 



modified-LP code, we adopt the interpolation scheme de- 
scribed in Appendix E 4 of Hamimeche & Lewis (2008), 
using a high value of the multipole Zmax to maximize the 
accuracy. This allows running the simulation several times 
without excessive consumption of CPU time and memory. 
We work under the null hypothesis that tensor modes are 
absent in the early Universe, so that the produced B-mode 
polarization is due only to the power transfer from the pri- 
mary scalar E- modes into the lens-induced B-modes. We 
choose Zmax ~ 6143 and interp-factor = 1.5, effectively the 
same resolution as HEALPix'' with pixelization parameter 
Nsidc = 2048, which corresponds to an angular resolution of 
1.72' (Gorski et al. 2005), with l2N'i^^ pixels in total. 
For comparison and testing, we also generate the cor- 
responding unmodified-LP and unlensed T, Q, U maps, by 
choosing the same positive seed in the LP parameter file 
(hereafter "unmodified-LP" stands for the results obtained 
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by lensing the sky with the original unmodified version of 
the LensPix code). 



3 SIMULATED LENSED MAPS 

The subtraction of the unlensed CMB maps from the cor- 
responding lensed CMB maps allows highlighting the dark- 
matter distribution which causes the gravitational defiec- 
tion of the CMB photons. In the upper and middle panels 
of Fig. 2 we show the resulting difference temperature and 
polarization maps for a particular seed choice and for the 
MS-modified-LP case only, since the arcminute-scale differ- 
ences with respect to the unmodified-LP case are not visible 
by the naked eye. It is worth noting how the distribution of 
the deflection-angle modulus (lower panel) reflects itself in 
the distribution of the temperature difference map (upper 
panel) and in the distribution of the map of the polarization 
modulus AP — ^/AQ^ + AU^ (middle panel), which has 
been obtained from the difference Q and U maps. Unfor- 
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Figure 4. Difference maps between the ips lensed and unlensed fields, obtained in the MS-modified-LP case (upper panel) and in the 
unmodified-LP case (lower panel). Units in mK. 



tunately, since the primary unlensed CMB is unknown, the 
difference maps cannot be directly observed. Anyway the 
simulated difference maps can help to physically and visu- 
ally understand how large scale correlations are imprinted 
on the CMB due to the large scale modes in the deflection 
field, and to catch effects that are not observable due to 
primary unlensed CMB. 

The temperature and polarization difference maps have 
the peculiarity of including all the information inferred from 
weak-lensing on the primary unlensed CMB. Nonetheless, 
it is well known that the lensed B-modes of polarization 
are more sensitive to the non-linear evolution of the cosmic 
structures than the T, Q and U modes separately (Lewis 
2005). The difference between the unmodified-LP and MS- 
modified-LP cases lies exactly at the non-linear level, since 
the lensing-potential realizations differ at multipoles / > 400 
in the spherical harmonic space for the two cases. More- 
over, from Fig. 1, we observe that the MS lensing poten- 
tial shows an excess of power on I > 2500 with respect to 
the CAMB approximation, and, as already pointed out, the 



non-linear dark matter evolution may enhance the level of 
non-Gaussianity present in the lensed CMB maps. 

Consequently, in order to verify if the different distri- 
butions of the lensing potential in the MS-modified-LP and 
unmodified-LP cases could visibly affect the lensed polar- 
ization distributions, we have constructed the full-sky maps 
of the scalar ^pE and pseudo-scalar potentials, which are 
related to the Q and U Stokes parameters as follows (Bunn 
et al. 2003): 



• + iU =66 {ipE + i'4^B) 



' — iU = 66 {'ijjE — i^ib) 



(2) 



(3) 



where the spin-raising 6 and spin-lowering 6 operators on 
the sphere are defined as (Newman & Penrose 1966) 



= — si 



(4) 
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Figure 5. Top panel: The dot-dashed light green hne represents 
the ratio between the lens-induced B power spectrum, which is 
obtained using the multipole range < lip < 6143 for the spher- 
ical harmonic coefficients of the lensing potential, and the spec- 
trum obtained using only < lip < 400. The solid blue line and 
the dashed red line represent these ratios in the T and E cases, 
respectively. Bottom panel: The same as the upper panel when we 
calculate these ratios for the temperature and polarization power 
spectra extracted from the difference T, Q and U maps, in the 
< ^ip < 6143 and < (jp < 400 cases, respectively. 



6 



(5) 



and s is the spin of the function to which the operator is 
apphed. The quantities tpE and tpB are directly related to 
the electric and magnetic types of polarization, since their 
spherical harmonic coefficients can be written in terms of the 
harmonic coefficients of the E- and B-modes, respectively 



^Ylr. 



-Y,[{l^2)\/{l + 2y]^'^aB,imYi„ 



(6) 



(7) 



The and ipB potentials are very useful for real space 
calculations and, exploiting the HEALPix routine SYN- 
FAST (Gorski et al. 2005), we have produced the corre- 
sponding maps as synthetic realizations, using the E, B 
spherical harmonic coefficients extracted (via the HEALPix 
routine ANAEAST) from the lensed Q, U simulated maps 



and multiplied by the prefactors of Eqs. (6)- (7), respectively. 
On high multipoles these prefactors have the asymptotic 
form [{I - 2)\/{l + 2)!]^/^ ~ consequently their effect 
is to suppress the small scale power with respect to the E- 
and B-moAe cases, and to make the large scale structure 
differences much more evident. 

In fact, looking at Figs. 3-4, we observe that the lensed 
ipE and ipB difference maps have a rfegree-scale distribu- 
tion which differs in the MS-modified-LP and unmodified- 
LP cases, even if the corresponding lensing-potential maps 
differ on arcminute scales (/ > 400). This result is more ex- 
pected in the _B-mode case since it is well-known that the 
non-linear density field introduces ~ 10% corrections to the 
lens-induced B-mode power on all scales (Lewis & Challi- 
nor 2006), so that multipoles I > 400 will affect the 7/)s-field 
realizations on larger scales too. The same is instead less 
obvious in the iJ-mode case. 

In order to understand the reason of the different large- 
scale distribution of these maps, we have also produced 
the lensed temperature, tps and tpB fields using a Gaus- 
sian lensing-potential distribution which includes only the 
scales that correspond to multipoles up to hp = 400 in the 
spherical harmonic domain (where the subscript "Ip" stands 
for lensing-potential). For comparison among the figures, we 
have used the same seed as for the fully (0 < hp < 6143) 
lensed maps, even if the results have been tested for differ- 
ent seeds. In Fig. 5 we show the ratios between the signals 
corresponding to the two cases, both for the lensed- and 
difference- maps. 

In particular, if we consider the lensed T- and i5-mode 
power spectra (top panel of Fig. 5) which include the contri- 
bution from the primary CMB, we do not observe any sub- 
stantial difference on the large scales between the two cases 
< hp < 400 and < hp < 6143, so that on large scales 
the lensed electric polarization appears not to be much af- 
fected by the non-linear scales in the same way as the lensed 
temperature is. The B case is different since we are working 
with the null hypothesis of vanishing primary magnetic po- 
larization, so that we are considering here only lens-induced 
B-modes. 

Nonetheless, if we consider the power spectra extracted 
from the difference T and i/)_e maps (where we are subtract- 
ing the primary CMB), and again take the ratios between 
the signals in the two cases < hp < 400 and < /jp < 6143 
(lower panel of Fig. 5), we see that on large scales the E- 
mode polarization, cleaned of the primary signal, gets power 
from the smaller scales. In particular we see that, in the E 
and B cases, multipoles larger than hp = 400 transfer about 
100% of the power to the low multipoles I < 100, i.e. on 
scales large enough for the ffat-sky approximation to be in- 
adequate. This transfer of power towards low multipoles is 
possible only for peculiar alignment angles between the un- 
lensed i?-mode and the lensing structures, i.e. for almost, 
but not quite aligned modes. This is analogous to the long- 
wavelength beat mode one obtains from superposing two 
oscillations with two almost equal frequencies. We do not 
observe the same effect when we consider the ratio between 
the signals extracted from the temperature difference maps 
in the two cases < hp < 400 and < /;p < 6143, where the 
transfer of power is less than 10% on / < 100. As we move to 
larger multipoles instead, the power transfer decreases and, 
for 400 < I < 2000, the lensed E trend starts converging 
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Figure 6. Top panel: Temperature power spectra {^K^) for the 
different cases described in the text. Bottom panel: Temperature 
power spectrum ratios with respect to the linear lensed case. 
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Figure 7. Top panel: E-mode power spectra {^K^) for the dif- 
ferent cases described in the text. Bottom panel: B-mode power 
spectrum ratios with respect to the hnear lensed case. 



to the lensed temperature one, while the power of the lens- 
induced B-modes goes on increasing with the multipoles, 
still fed by the non-linear scales in the lensing-potential at 
lip > 400. 

Therefore we conclude that in Figs. 3-4 the differ- 
ent degree-scale distributions which characterize the MS- 
modified-LP and unmodified-LP t/jb {'4'b) difference map re- 
alizations (even though the corresponding lensing-potential 
realizations differ only on Up > 400) represent the imprints 
of the power transfer from smaller to larger scales induced 
by the lensing remapping onto the lensed E and B signals 
cleaned of the primary CMB, in a way which differs from 
what happens to the lensed temperature. 



4 ANGULAR POWER SPECTRA 

In this Section we perform several quantitative analyses of 
the obtained results. Our principal aim is to test the consis- 
tency and quantify the differences between the theoretical 
expectations and the findings from our simulation proce- 
dure. 

As a first check, we have extracted from the cor- 
responding maps the angular power spectra of the MS- 
modified-LP and unmodified-LP simulated lensed T-, E- 



and B-components, together with the TE lensed cross- 
correlation power spectrum. This is simply done with the 
use of ANAFAST, adopting the correct deconvolution rules 
for the HEALPix pixel window functions, and with Eqs. (6)- 
(7) relating the spherical harmonic coefficients of the ipE and 
tpB potentials to the corresponding coefficients of the E- and 
B-modes. 

We have also checked that the MS-modified-LP and 
unmodified-LP simulated spectra extracted directly from 
the maps agree with the ones obtained as direct outputs 
of the MS-modified-LP and unmodified-LP codes, respec- 
tively. Moreover, to test the accuracy, we have compared all 
the angular power spectra with the lensed CMB power spec- 
tra obtained using the all-sky correlation function technique 
implemented in CAMB (Challinor & Lewis 2005). These an- 
gular power spectra are shown in the upper panels of Figs. 6- 
9 respectively, together with the CAMB unlensed and the 
linearly lensed spectra, where the latter are obtained in the 
linear approximation of the cosmic structure growth. More- 
over, in the lower panels of the same figures we show the 
ratios between the simulated lensed CMB angular spectra 
with respect to the corresponding linear signal. 

In the TT, EE and TE cases, it is clearly visible that 
CMB lensing smears out the acoustic peaks by transferring 
power from larger to smaller scales. Moreover, since we are 
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Figure 8. Top panel: Lens-induced B-mode power spectra {^K^) 
for the different cases described in the text. Bottom panel: Lens- 
induced B-mode power spectrum ratios with respect to the linear 
lensed case. 



working under the null hypothesis of vanishing intrinsic ten- 
sor modes, Fig. 8 represents the power spectrum of the lens- 
induced B-modes into which part of the primary i5-modes 
has been converted as a result of the displacements and dis- 
tortions induced from the gravitational deflection onto the 
electric-type polarization field. 

As Figs. 6, 7 and 9 show, for all the lensed TT, EE 
and TE angular spectra, and for all multipole orders up to 
I < 3500, we observe a mostly perfect agreement between the 
MS-modified-LP simulated signals and the unmodified-LP 
and CAMB ones, where the non-linear structure evolution 
is semi-analytically taken into account. On the very large 
multipoles I > 3000, a small difference appears between the 
CAMB and unmodified-LP lensed TT, EE and TE spec- 
tra, probably due to a numerical effect deriving from the 
different computational machinery implemented in CAMB 
and LensPix respectively. The non-linear effects start to be 
important at / > 2500, according to the semi-analytical ex- 
pectations, and, on multipoles I ~ 4000, which correspond 
to angular scales of few arcminutes, they grow up to ~20% 
for the lensed temperature and up to ~10% for the lensed 
electric polarization. Moreover, the excess of power present 
in the MS lensing-potential power spectrum at multipoles 
I > 2500 (see Fig. 1) manifests itself as a slight excess of 
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Figure 9. Top panel: T-E cross-power spectra {fiK^) for the 
different cases described in the text. Bottom panel: The fractional 
change in the lensed T-E cross-power spectrum due to the non- 
linear matter evolution for the different cases described in the 
text. 



power in the MS-modified-LP lensed TT and EE spectra 
at Z > 3500. In particular, at I — 4096, the non-linearities 
present in the Millennium Simulation produce a ~6% ex- 
cess in the lensed temperature and a ~1.3% excess in the 
lensed electric polarization, with respect to the unmodified- 
LP case. 

On the other hand, the excess of power due to the MS 
non-linearities is much more evident in the angular power 
spectrum of the magnetic component (see Fig. 8) where this 
effect starts already at I ^ 2500 and grows up to ~14% at 
I = 4096, with respect to the unmodified-LP case. Moreover, 
the non-linear effects, which in the lensing potential appear 
at I > 400, spread on all the scales of the lens-induced B- 
mode spectrum, being already of the order of ~7.5% on 
all the multipoles I < 1000, and growing up to more than 
~70% at I = 4096. As we have previously noticed, this non- 
linear effect on the magnetic polarization is simply explained 
considering that the lens-induced B-modes are very sensitive 
to the non-linear evolution of the cosmic structures. 

For what concerns the non-linear effects, as noticed al- 
ready in Sec. 3, it is important to stress again that, when 
considering the power spectra extracted from the difference 
maps (in which the primary CMB has been subtracted) , we 
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Figure 10. Top panel: The dot-dashed Ught green Une represents 
the ratio between the lens-induced B power spectrum, which is 
obtained including the non-Unear cosmic structure evolution, and 
the spectrum obtained in the linear limit. The solid blue line and 
the dashed red line represent these ratios in the T and E cases, 
respectively. Bottom panel: The same as the upper panel when we 
calculate these ratios for the temperature and polarization power 
spectra extracted from the difference T, Q and U maps, in the 
linear and non-linear cases, respectively. 



notice a very similar large-scale behavior between the elec- 
tric and magnetic polarization, which strongly differs from 
the temperature trend. In order to analyse this effect, we 
have produced the /meari?/-lensed T, Q, U maps and sub- 
tracted the unlensed field. In Fig. 10 we show the two cases: 
the ratios of the lensed T, E, B spectra with respect to 
the corresponding spectra extracted from the linearly-lensed 
maps, and the same ratios extracted from the corresponding 
lensed and linearly-lensed difference maps. While in the first 
case (top panel) the transfer of power from small scales to 
large scales is only visible in the lens-induced B-mode spec- 
trum, in the second case (bottom panel) this power transfer 
is also observable in the lensed iJ-mode spectrum, with the 
same amount of ~7.5% up to i < 100 as for the lens-induced 
_B-modes. The same does not occur to the lensed tempera- 
ture even after subtracting the unlensed CMB. Finally, mov- 
ing to higher multipoles, the E and T signals converge to 
the same trend, while the B signal keeps on growing because 
of the non-linear power. 



5 ONE-POINT STATISTICS 

As a second analysis, we compare the one-point probability 
density distributions (PDFs) of the lensed T, Q, U maps 
with the PDFs of Gaussian distributions randomly gener- 
ated using the corresponding mean and standard deviation 
values of the lensed maps. The PDFs of both the simu- 
lated MS-modified-LP and unmodified-LP T, Q, U maps do 
not show evidence of non-Gaussian features (upper panel of 
Fig. 11), and are consistent with the unlensed temperature 
and polarization distributions within a few per cent. This oc- 
curs thanks to the high angular resolution of our maps (see 
e.g. Kesden et al. (2002)), and is understandable because 
the main effect of CMB lensing is to transfer power among 
different scales, without generating new power (and we as- 
sume that the primary unlensed CMB field has a Gaussian 
distribution) . 

On the other hand, it is well-known that the remap- 
ping induced by lensing onto the T, Q, U fields generates a 
non-Gaussian signature in the lensed sky which is optimally 
characterized by higher order statistics as the bispectrum 
and trispectrum (see e.g. Lewis & Challinor 2006, and ref- 
erences therein). This is mainly due to the fact that the 
lensed CMB can be considered as a function of two fields, 
the unlensed sky and the lensing potential, which in first ap- 
proximation can be assumed to be Gaussian (unmodified-LP 
case), even if the non-linear evolution of large-scale struc- 
tures produces non-Gaussian features in the distribution 
of the projected potential (MS-modified-LP case). In par- 
ticular, the non-Gaussianity associated with the large-scale 
structure induces non- Gaussian contributions to the distri- 
bution of the lensing potential such that its n-point correla- 
tor in Fourier space will be non- vanishing for some value of 
n, and this, in turn, has an impact on the connected part of 
the n-point correlators of the lensed CMB (see e.g. Kesden 
et al. (2002)). Anyway, there exist also other phenomena 
that can produce non-Gaussian effects in the lensed CMB, 
for instance the correlations between CMB lensing and the 
Sunyaev-Zeldovich effect. 

Here, as a first characterization of the non-Gaussianity 
strength produced by CMB lensing, we consider the PDFs 
of the difference T, Q, U maps. In this way, we subtract the 
unlensed Gaussian sky and isolate the non-Gaussian term 
which physically generates the non-Gaussian signatures in 
the lensed CMB. The obtained PDF in the temperature case 
is showed in the lower panel of Fig. 11 (the Q and U cases 
have a similar trend). 

The PDFs of the difference maps are characterized by a 
kurtosis excess with respect to Gaussian distributions ran- 
domly generated with the same mean and standard devi- 
ation values. More precisely, averaging over different real- 
izations, the kurtosis excess is ~ 2.41 for the difference T- 
map, ~ 2.14 for the difference Q-map and ~ 2.16 for the 
difference [/-map. Actually we would expect that the ex- 
cess of non-linearities in the dark-matter distribution of the 
Millennium Simulation should show up as an excess of non- 
Gaussianity in the MS-modified-LP lensed difference T, Q, 
U maps, with respect to the unmodified-LP case. However, 
the one-point statistic is probably unable to capture this ef- 
fect, and Fig. 11 shows in fact that there is not a significant 
difference between the unmodified-LP and MS-modified-LP 
cases. In this respect, we believe that a more suitable es- 
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Figure 11. Top panel: The PDF of the temperature maps ob- 
tained in the MS-modified-LP and unmodified-LP cases compared 
to the PDF of Gaussian distributions with the corresponding 
mean value and standard deviation, respectively. Bottom panel: 
The same as in the upper panel for the PDFs of temperature 
difference maps. 



timator should be developed in order to detect the contri- 
bution from the non-linear structure evolution to the total 
non-Gaussian statistics of the lensed maps. We reserve this 
analysis for future work. 



6 CONCLUSIONS 

We have constructed the first all-sky CMB temperature and 
polarization lensed maps based on a high-resolution cosmo- 
logical A^-body simulation, the Millennium Simulation. 

To this purpose we have exploited the lensing potential 
map obtained using the map-making procedure developed 
in Carbone et al. (2008) which integrates along the line-of- 
sight the MS dark matter structures by stacking and ran- 
domizing the simulation boxes up to z = 127. Specifically, 
we have modified the LensPix code (Lewis 2005) by sup- 
plying it with the spherical harmonic coefficients extracted 
from the MS lensing potential map and by implementing 
directly in the code itself the large-scale structure adding 



technique (see Carbone et al. 2008) which allows to rein- 
state the large-scale power in the angular lensing potential 
that is not correctly sampled by the A^-body simulation. In 
this way we also preserve the correct correlation between the 
lensing potential and the ISW effect on multipoles I < 400 
in the simulated temperature anisotropics. 

Using our modified version of the LensPix code, we have 
constructed lensed CMB simulated maps with ~ 5 million 
pixels and an angular resolution of ~ L72', based on poten- 
tial fields calculated on 2560"^ mesh cells from the Millen- 
nium Simulation. 

After subtraction of the unlensed maps, the correspond- 
ing lensed T, Q, U simulated maps refiect clearly the same 
large scale structure which is present in the modulus of the 
angular gradient of the lensing potential map (see Fig. 2). 

We have also constructed the maps of the scalar and 
pseudo-scalar potentials tpE and ^ps which are directly re- 
lated to the electric and magnetic types of polarization, re- 
spectively. Their difference maps, shown in Figs. 3-4, present 
distinct degree-scale distributions in the MS-modified-LP 
and unmodified-LP cases (where the latter is obtained from 
the unmodified LensPix code as previously explained), ow- 
ing to the power transfer from non-linear to large scales in 
the lensed E and B fields. 

As a quantitative study of the simulated maps, we have 
performed power spectrum and one-point statistics analyses. 
We find that the lensed TT, TE, EE and BB power spec- 
tra, obtained using the MS dark matter distribution, mostly 
overlap with the corresponding semi-analytic expectations 
on a range of multipoles up to Z ~ 2500 (see Figs. 6-8). 
This latter result points out the effectiveness of the Mil- 
lennium Simulation in reproducing correctly the findings of 
the theoretical approach. Furthermore, an excess of power 
is observable in the MS-modified-LP case on larger multi- 
poles, in particular for the BB spectrum. We believe that 
this excess originates from the accurate inclusion of non- 
linear power in the Millennium Simulation, which is present 
in the MS lensing potential power spectrum as well (Fig. 1). 
This outcome should be taken into account in the various 
delensing approaches, since these non-linear effects can have 
an impact on the quality of the reconstructed unlensed sky, 
particularly in view of the detection of primary B-modes. 

Finally, we have derived one-point statistics both of the 
simulated T, Q, U lensed maps and of the corresponding 
difference maps, in the MS-modified-LP and unmodified-LP 
cases. The comparison of the resulting PDFs with respect 
to Gaussian random distributions with the same mean and 
standard deviation does not show any statistical difference 
when the total lensed maps are considered. On the other 
hand, we find that the difference maps are characterized 
by a kurtosis excess. This result represents the distinctive 
weak-lensing effect of inducing non-Gaussianity in the un- 
lensed Gaussian CMB field. This is simply explained if in 
first approximation we think of the lensed field as the prod- 
uct of two Gaussian fields, i.e. the lensing potential and 
the primary unlensed CMB. Actually, correctly speaking, 
the lensing potential map, derived by integration of the MS 
dark matter distribution, preserves an intrinsic degree of 
non-Gaussianity due to the non-linear evolution of the cos- 
mic structures, which at some level contributes to the non- 
Gaussian statistics of the simulated lensed CMB. Unfortu- 
nately the one-point statistic is not sufficient to disentangle 
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the two contributions, i.e. the dominant effect coming from 
the product of two different fields, and the subdominant ef- 
fect coming from the non-linear matter evolution. This is 
the reason why in Fig. 11 there seems to be no difference 
between the MS-modified-LP and unmodified-LP cases. We 
reserve a more accurate characterization of non-Gaussianity 
in CMB lensing statistics to a future work. 
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